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Abstract 

Incentives, problems and solutions associated with 
interfacing laboratory equipment with multiuser, virtual 
memory computers are presented. The major difficulty 
concerns how to utilize these computers effectively in a 
medium sized research group. This entails optimization of 
hardware interconnections and software to facilitate multiple 
instrument control, data acquisition and processing. The 
architecture of the system that was devised and associated 
programming and subroutines are described. An example 
program involving computer controlled hardware for 
ultrasonic scan imaging is provided to illustrate the operational 
features. 


Introduction 

Computers based on “virtual memory”, multiuser 
architecture have recently become economically available to 
medium sized research groups. These new computers offer 
advantages to the experimentalist, though acceptance has been 
sometimes limited. One of the most important advantages of 
a “virtual memory” machine is the capability of fast access 
to arrays having dimensions greater than 64 kbytes. Arrays 
larger than 64 kbytes were available on some older computers 
based on “physical memory.” The access time to the memory 
beyond 64 kbytes is extremely slow for these systems. Fast 
access to large arrays is particularly important now where large 
data sets are required for research systems. 

Many researchers are accustomed to a dedicated computer 
running a single-user operating and data acquisition system. 
When there are several experimentalists in a research group, 
the advantage of a multiuser system becomes important. The 
resources a computer provides can be shared among the users 
at a much lower cost than that of maintaining a separate 
dedicated computer and peripherals for each project. Software 
routines for instrument control and data acquisition as well 
as data bases are sharable resources. Expensive hardware 
resources, such as array processors, images processors, 
terminals, printers and color plotters can also be shared. 

There are difficulties in assembling the hardware and 
software facilities required for a multiuser, experimental 
research oriented system. The single-user operating systems 
typically use built-in commands for easy instrument 


communication and other functions well suited for the 
experimenter’s needs. A multiuser computer’s input-output 
capabilities are more powerful but harder to understand. The 
user is left with the responsibility for bridging the gap between 
the capabilities provided and his particular needs for data 
acquisition and communication with devices. 

The general purpose interface bus (GPIB) is used extensively 
in single-user computer systems for data acquisition and 
instrument control. For single-user systems the software for 
data acquisition and instrument control is often provided by 
the manufacturer of the system. This same bus can be used 
in a multiuser system. However, the necessary software is not 
provided. This lack of software has been a major obstacle for 
researchers and has prevented them from moving from a single 
to a multiuser virtual memory system. 

In this work we will provide the reader with the information 
required for developing a multiuser computer system that 
meets the needs of a medium sized research group. The basic 
computer system and its components are described. Specific 
examples of GPIB instrument control and data acquisition with 
a multiuser system are described for use in both BASIC and 
FORTRAN programming. Additionally, a collection of high 
and low level subroutines for computer control of instruments 
have been listed in the appendixes. 

Computer System and Instruments 

Computer System 

A computer system for a research group has several 
components. Figure 1 shows a block diagram of a balanced 
computer system specifically designed for an experimentally 
intensive research group. The system is based on a virtual 
memory computer. There are three central processing units 
(CPUs) in this configuration, each with 16 Mbyte memory, 
and two, 474 Mbyte, hard disks attached. The second hard 
disk is for backup of system and user software. Generally 16 
Mbytes of memory is required for the processing of large data 
arrays. Each CPU (subsystem) has a major specific function. 
CPU 1 is designed for data analysis and contains an array 
processor for fast Fourier transformations. In addition, CPU 
1 shares an image processor with CPU 2 for display images 
viewed with the image processor. CPU 2 is designed for image 
analysis and display. The 32 Gbyte optical disk drive provides 
easy access and storage of images. CPU 3 is designed for data 
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Figure 1.— Computer system. 


acquisition and contains four GPIB interfaces yielding eight 
GPIB channels (ports) and an analog to digital (A/D) 
input/output interface with 16 channels. Up to nine users can 
use the third CPU simultaneously for data acquisition. Of 
course A/D interfaces may be added for additional users. 
Alternatively, fewer users could run multiple experiments 
simultaneously. 

It is important to note that, although each CPU subsystem 
has been designed for a specific function, it is not limited to 
that function. Therefore, CPUs 1 and 2 also have GPIB and 
A/D interfaces. As shown a total of 15 researchers may acquire 
data or control instruments simultaneously via the GPIB or 
A/D interfaces. Equally important to the number of users 
working with the instruments is that additional users may be 
simultaneously analyzing data, image processing, developing 
software and generating reports. 

The three CPUs are connected via an Ethernet connection 
for easy transfer of files between systems. A terminal server 
connected to the Ethernet allows for multiuser, remote access 
to each of the CPUs. This remote access allows researchers 
to have terminals and instruments placed at different locations 
throughout the site. Off-site access to the system is permitted 
via modems attached directly to the CPU or Ethernet. Modems 
are used for initiating, stopping or checking the status of an 
experiment on evenings or weekends. 


A variety of instruments may be attached to the GPIB 
interface. It is this interface that is often the stumbling block 
in meeting the researchers’ needs. The problem is two-fold 
for a multiuser system. First, the GPIB bus has a limit on total 
length of 20 m. One cannot expect to have all experiments 
taking place in one location (i.e., room, building, etc.); 
therefore, the need arises to extend the GPIB bus beyond this 
limit. Second, friendly software for multiuser, virtual memory, 
GPIB based instrument control and data acquisition is not 
available. 

A typical configuration for a nondestructive evaluation 
research group is shown in figure 2. Several research 
instruments including ultrasonic, x-ray tomographic, and 
tunneling scanners are attached to the computer system. 
Recently, GPIB bus extenders have been made available from 
several manufacturers. By using bus extenders, these 
instruments may be placed at locations remote from the CPUs. 
The GPIB bus extenders are transparent to the researchers’ 
software so that no software modifications are needed. 

The software required to control these instruments and for 
data acquisition is the most complicated of the two problems 
and for this reason it will be presented in detail. 

Programming the components of the precision acoustic 
scanner covers a wide range of software problems encountered 
when connecting instruments to a computer via the GPIB bus. 
We will focus on this system, as it represents the most common 
problems. 



Figure 2.— Instrument configuration for a nondestructive evaluation 
research group. 
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The components of the precision acoustic scanning system 
are shown in figure 3. The specimen to be interrogated is 
mounted on a X-Y-Z positioning table. In contact with the 
sample is a ultrasonic transducer mounted on a displacement 
pressure gauge. During operation, the positioning table is 
moved to a series of points that describe an array. This is done 
while maintaining a constant pressure between the transducer 
and specimen. At each point ultrasonic waveforms are 
selected, digitized and stored for later analysis. 

This system has the following six components that 
communicate over the GPIB bus: three positioning tables, 
waveform digitizer with a time base and voltage base, time 
delay and voltmeter. The programming for this rather complex 
scanning system contains all of the features needed for 
adaptation to other computer controlled systems. 



X-Y-Z POSITIONING 
TABLE 


WAVEFORM DIGITIZER 


TIME BASE 
VOLTAGE BASE 


PRESSURE SENSOR 
VOLTMETER 


TIME DELAY 


ULTRASONIC 

PULSER- 

RECEIVER 




Figure 3.— Precision acoustic scanning system. 


Progr ammin g for Instrument Control and 
Data Acquisition 

In this section a fundamental program, written in 
FORTRAN, is presented that operates the precision acoustic 
scanning system. The subroutines for this program are 
provided in the appendixes. 

Initially a GPIB channel is assigned to the precision acoustic 
scanner. Here a channel represents a physical GPIB port on 
the CPU. 

The FORTRAN statement 


CALL STRTGPIB (0) 

initializes channel 0. Any instrument addressed subsequently 
in the program will communicate through channel 0. Next each 
component of the scanning system must be initialized. There 
are six components, and their physical addresses have been 
adjusted to be contiguous from 32 to 37. These are also called 
addresses 0 to 5 where the it is implied that the standard base 
GPIB address of 32 is to be added. The physical addresses 
for the digitizer, time delay, X-axis, Y-axis, Z-axis and 
voltmeter are set by the researcher at 0, 1, 2, 3, 4 and 5, 
respectively. 

DO 10 NUMBER = 0, 5 
CALL INITINSTR (NUMBER) 

10 CONTINUE 


As indicated earlier, the specimen is scanned in an organized 
array in the X-Y plane. The step size in the X, Y and Z 
directions are stipulated by the user as 


XSTEP = 100 
YSTEP = 100 
ZSTEP = 10 


the z axis will not be used for this 
template. But an argument is 
required for these subroutines. 


and sent to the X-, Y- and Z-axis of the positioning table 


CALL SETXYZ (2, XSTEP, 3, YSTEP, 4, ZSTEP) 

where 2, 3 and 4 correspond to the X, Y and Z positioner 
axes. The ultrasonic waveform needs to be digitized and 
stored. The ultrasonic wave occurs at specific time and is set 
by the time delay subroutine 

CALL SETDELAY (1, DELAY(l)) 
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where the first argument is the address of the time delay. 
DELAY(l) is the time delay required to place the ultrasonic 
wave within the digitizing window. The volts and time per 
division of the time and voltage base are set 

CALL PUTTIME (0, SETTING) 

CALL SETVOLTDIV (0, SETTINGS(l)) 

where the first argument is the primary address of digitizer 
that contains the time and voltage bases. The SETTING and 
SETTINGS(l) are the proper time and voltages per division 
for displaying the ultrasonic wave in the digitizing window. 
The actual digitization of the ultrasonic waveform occurs with 
the command 

CALL GETSA (0, 64, 320, A) 

where the first argument is the address of the digitizer and 
64 is the number of samples that are to be averaged, 320 is 
the intensity of the writing beam in the digitizer and “A” is 
the array that will contain the digitized waveform. The 
waveform is saved to a file opened as unit 6 with 

WRITE (6) A 

After a waveform is collected at one point, another point 
is selected with the X-Y-Z positioning table with 

CALL MOVXYZ (2, 1, 1) 

where the first argument corresponds to the address of the X- 
axis. The second argument corresponds to the axis to be moved 
one step (1,2 and 3 for the X-, Y- and Z-axis, respectively). 
If the third argument is greater or less than zero, then the 
motion is in the positive or negative direction, respectively. 

This program would loop back to the CALL PRESSURE 
statement to complete a scan line. The above FORTRAN 
source code gives a template for writing other source code. 
The complete source code for scanning a 5- by 5-mm square 
area is given in appendix A. 

Appendix B contains a listing and description of each of the 
user callable high level subroutines used for instrument control. 
The low level source code for these routines is provided in 
appendix C. This source code may be adapted for writing other 
instrument control subroutines for additional instruments not 
covered in this work. 


The complete program SCAN. FOR must be compiled first 
$ FORTRAN SCAN 
and linked 

$ LINK SCAN, SYSSSYSROOT : [SYSLIB]IMAGELIB/LIB - 
SYS$SYSROOT:[SYSLIB]IEXSUI.OBJ 

then run with 

$ RUN SCAN 

The subroutines provided as well as this and other data 
acquisition and instrument control programs may be shared 
with other users. Researchers may also write programs for 
operating instruments attached to other GPIB channels. With 
this capability, several researchers may be operating 
instruments and acquiring data simultaneously. A considerable 
amount of time is saved by sharing these developed programs. 

We have found that, when programs and instruments are 
easily shared, program development occurs at a rapid rate. 
This rapid growth occurs because researchers modify each 
control program to meet their particular needs. In doing so, 
the fundamental program expands and is fine tuned by all the 
researchers to form a very versatile system. 


Summary 

A multiuser, virtual memory computer system was devised 
and programmed for use by a medium sized research group. 
Problems associated with optimizing the system for shared, 
multiuser applications were resolved. Particular attention was 
given to hardware interfacing and development and 
compilation of subroutines. The result is a straightforward 
operating network for instrument control, data acquisition and 
processing. A comprehensive example program is provided 
for illustrative purposes. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 17, 1989 
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Appendix A 

FORTRAN Program for Running Precision Acoustic Scanning System 


C PROGRAM SCAN. FOR (November 7, 1988) 

C by Edward R. Generazio and David B. Stang 

C 

C WRITTEN IN DEC FORTRAN VERSION, COMPILED AND RUN ON A DEC MICROVAX II 
C WITH MVMS VERSION 4.5 OPERATING SYSTEM 
C 

C THIS PROGRAM RUNS THE PRECISION ACOUSTIC SCANNING SYSTEM DEVELOPED AT NASA 
C LEWIS RESEARCH CENTER, STRUCTURAL INTEGRITY BRANCH, STRUCTURES DIVISION, 

C CLEVELAND, OHIO. 

C 

C THE INSTRUMENTS ATTACHED ARE 
C 

C TWO KLINGER SCIENTIFIC C-1.22 PROGRAMMABLE STEPPER 

C MOTOR CONTROLLER 

C 

C TEKTRONIX 7912 PROGRAMMABLE DIGITIZER WITH A TEKTRONIX 

C 7B90P PROGRAMMABLE TIME BASE AND 

C A TEKTRONIX 7A16P PROGRAMMABLE AMPLIFLIER (VOLTAGE BASE) 

C 

C HEWLETT-PACKARD 5359A TIME SYNTHESIZER (TIME DELAY) 

C 

C FLUKE 8520A DIGITAL MULTIMETER (VOLTMETER) 

C 

C SENSOTEC LOAD SYSTEM. THE SENSOTEC ANALOG OUTPUT IS PROPORTIONAL 

C TO THE PRESSURE ON THE LOAD CELL AND IS MEASURED BY THE 

C VOLTMETER. 

C 

C 

C 

C THIS PROGRAM MOVES THE SCANNER IN AN ORGANIZED 5MM BY 5MM 

C SQUARE ARRAY . 

C 

C AT EACH POINT THE CONTACT PRESSURE BETWEEN THE TRANSDUCER 

C AND THE SPECIMEN BEING SCANNED IS 15 TO 17 POUNDS. 

C 

C ALSO AT EACH POINT TWO WAVEFORMS ARE DIGITIZED AND STORED IN 

C MASS STORAGE FOR LATER ANALYSIS. 

IMPLICIT NONE 

INTEGER*4 XSTEP , YSTEP , ZSTEP , XAXIS , YAXI S ! DEFINITIONS 
INTEGER* 4 I , J , K ,NAVE , MAI 
INTEGER*2 A(512) 

BYTE BUFFER ( 32 ) ,WRKO( 80 ) 

REAL DELAY ( 2 ) , VOLTSET( 2 ) 

CHARACTER* 10 FILNAME 

INSTRUMENTS ARE SET AT FOLLOWING ADDRESSES 

DIGITIZER - 0 
TIME DELAY - 1 
X AXIS - 2 

Y AXIS - 3 

Z AXIS - 4 

VOLTMETER - 5 
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FILNAME- ' SCANNER . DAT ' ! NAME OF FILE TO SAVE WAVEFORM IN 

OPEN ( UNI T-6, STATUS-' NEW' , NAME- ' FI LNAME ' , FORM- ' UNFORMATTED ' ) 

I OPEN THE FILE 

DELAY ( 1 ) «4 . 3 5E-6 ! TIME DELAY FOR FIRST WAVEFORM 

DELAY ( 2 ) «*5 . 20E-6 ! TIME DELAY FOR SECOND WAVEFORM 

VOLTSET( 1 )=0 . 5 ! VOLTS PER DIVISION FOR FIRST WAVEFORM 

VOLTS ET( 2 )=0 . 2 ! VOLTS PER DIVISION FOR SECOND WAVEFORM 

XSTEP-1000 ! STEP SIZE FOR THE X AXIS 

YSTEP-1000 ! STEP SIZE FOR THE Y AXIS 

ZSTEP-10 1 STEP SIZE FOR THE Z AXIS 

XAXIS-5 ! NUMBER OF STEPS ALONG THE X AXIS 

YAXIS-5 ! NUMBER OF STEPS ALONG THE Y AXIS 

NAVE-64 ! NUMBER OF SAMPLES TO AVERAGE 


CALL STRTGPIB ( 0 ) 1 INITIALIZE CHANNEL 0 OF THE GPIB 

! INTERFACE 

DO 10 1=0,5 I INITIALIZE INSTRUMENTS AT GPIB 

10 CALL INI TINSTR ( I ) ! ADDRESSES 0,1, 2, 3, 4 AND 5 

CALL SETXYZ ( 2 , XSTEP , 3 , YSTEP , 4 , ZSTEP ) 

! TELL THE POSITIONER THE SIZE OF THE 
! STEPS TO TAKE 

CALL GETMAI ( 0 , MAI ) ! GET THE CURRENT INTENSITY SETTING 

! ON THE DIGITIZER 

DO 20 I— 1 , YAXI S ! THE YAXIS WILL MOVE 5 TIMES 

DO 30 J— 1 , XAXI S ! THE XAXIS WILL MOVE 5 TIMES 

CALL PRESSURE ! ADJUST THE PRESSURE TO THE 

! APPROPRIATE VALUE 

DO 40 K-1,2 ! DO THIS FOR EACH WAVEFORM 

CALL PUTTIME ( 1 , DELAY ( K ) ) ! SET THE DELAY TIME FOR THE WAVEFORM 

CALL SETVOLTDI V ( 0 , VOLTS ET ( K ) ) 

! SET THE VOLTS PER DIVISION FOR 
! THE WAVEFORM 

CALL GETSA ( 0 , NAVE , MAI , A ) 

! SET THE DIGITIZER INTENSITY TO MAI 
l AVERAGE OVER 64 SAMPLES 
! PUT DIGITIZED WAVEFORM IN ARRAY 
! CALLED "A" 

40 WRITE ( 6 ) A ! SAVE WAVEFORM IN FILE NAMED SCANNER . DAT 

IF (J.EQ. XAXIS) GOTO 30 ! IF THE END OF A SCAN LINE IS REACHED 

! RETURN TO THE X AXIS ORIGIN AND 
! INCREMENT THE Y AXIS 
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30 


20 


CALL MOVXYZ (2,1,1) 
CONTINUE 

IF(I.EQ.YAXIS) GOTO 20 

CALL MOVORG ( 1 , XAX I S ) 
CALL MOVXYZ (3, 2,1) 
CONTINUE 
CLOSE ( 6 ) 

END 


! OTHERWISE MOVE THE X AXIS ANOTHER STEP 


! THE SCAN IS OVER END THE PROGRAM 
! CLOSE THE SCANNER . DAT FILE AND 
! RELEASE THE GPIB CHANNEL 

! OTHERWISE RETURN TO THE X AXIS ORIGIN 
! AND 

! MOVE THE Y AXIS ANOTHER STEP 


! CLOSE THE FILE SCANNER . DAT 
! END THE PROGRAM SCAN. FOR 


INCLUDE ' GPIBASE . FOR' 


1 INCLUDE THE INSTRUMENT CONTROL 
! AND DATA ACQUISITION SUBROUTINES 
! IN THE FILE BASE. FOR 
! WHEN COMPILING. 

1 THESE SUBROUTINES ARE IN THE 
! APPENDIX OF THIS ARTICLE 


TWO SUBROUTINES THAT CALL SUBROUTINES IN BASE. FOR 


SUBROUTINE PRESSURE 


ADJUST 

<<< 

Z AXIS TO 
+1 = DOWN 1 

GET A 
>>> 

GOOD PRESSURE 


PUPPER » 

.85 

j 

UPPER PRESSURE = 17 

lbs 

PLOWER - 

.75 

j 

LOWER PRESSURE =15 

LBS 

CALL GETFLUKE ( 5,P 

) | 

GET THE OUTPUT VOLTAGE FROM THE 
PRESSURE SENSOR 


I F ( P.LT.PLOWER ) CALL MOVXYZ ( 4 , 3, 1 ) 

! IF PRESSURE IS TO LOW THEN INCREASE 
! THE PRESSURE BETWEEN THE SPECIMEN 
! AND THE TRANSDUCER 


IF( P.GT.PUPPER ) CALL MOVXYZ(4, 3, -1 ) 

! IF PRESSURE IS TO HIGH THEN DECREASE 
1 THE PRESSURE BETWEEN THE SPECIMEN 
I AND THE TRANSDUCER 


IF( P.GE.PLOWER .AND. P.LE.PUPPER )GOTO 900 



! THE PRESSURE IS WITHIN AN 

ACCEPTABLE 


! RANGE 


GOTO 200 

• CHECK PRESSURE 


RETURN 

l PRESSURE ADJUSTED- RETURN 

TO MAIN PROGRAM 

end 
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SUBROUTINE MOVORG< SCANAXIS, NSTEPS ) 

MOVE SCANNER BACK TO THE SCANAXIS ORIGIN 

INTEGERS SCANAXIS 
INTEGER*2 NSTEPS 


DO 1000 K=l,1000 
1000 CALL MOVXYZ ( 4 , 3,-1 ) 

DO 1200 K=1 , NSTEPS-1 
1200 CALL MOVXYZ ( 2, SCANAXIS, -1 ) 
DO 1400 K-1,1000 
1400 CALL MOVXYZ ( 4, 3, 1 ) 

CALL PRESSURE 

CALL MOVXYZ ( 2, SCANAXIS, -1 ) 
CALL MOVXYZ ( 2, SCANAXIS, 1 ) 
RETURN 
end 


! MOVE Z AXIS UP 

! MOVE X AXIS BACK 

! MOVE Z AXIS DOWN 
! ADJUST PRESSURE 
I MOVE BACK AND FORTH 
! TO MAKE A GOOD CONTACT 
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Appendix B 

User Callable High Level Subroutines for Instrument Control 


GPIBASE subroutines 


GPIBASE refers to a DEC VAX/VMS object code set containing 
FORTRAN callable subroutines which can be used to communicate with 
various GPIB ( General Purpose Interface Bus, also known as 
IEEE-488 ) instruments found in the Structural Integrity Branch 
laboratories. The "DEC IEX Interface" consists an hardware 
driver, version 3.0, along with software which is documented 
in the "IEX-VMS-DRIVER" manual. The following describes 
subroutines which utilize the IEX Interface to perform 
tasks typically done in the labs. 

Only one IEX channel is operated per main 
program, and each subroutine requires an integer argument 
which is the GPIB address of the instrument to be 
addressed. For instance, the following is a table of the 
configuration in the Stuctural Integrity Laboratories: 


CHANNEL, GPIB ADDRESS INSTRUMENT 


0 , 0 
0 , 1 

0 , 2 

0, 3 
0, 4 
0, 5 


Tektronix 7912AD Digitizer 
Hewlett Packard 5359A Time 
Synthesizer 1 
Klinger Scientific c.122 
XYZ Controller- X Axis 
Y Axis 
Z Axis 

Fluke 8520A Multimeter 


1 , 1 
1 , 2 

1, 3 

1, 5 


Tektronix 7854 Digital Oscilloscope 
Testech LS86-C Ultrasonic Immersion 
Scanner 

Hewlett Packard Time 
Synthesizer 2 
Hewlett Packard 53591A 
A/D Converter 


GENERAL GPIBASE SUBROUTINES 


STRTGPIB(N) - Assign and initialize IEX channel N 

This must be called before calling channel N instrument 
subroutines. N is less than M and M is the number of 
physical ports. It assigns the VAX/VMS internal address 
and also clears the interface. Any instrument addressed 
subsequently in the program will communicate through 
channel N. 
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INITINSTR( ADDR ) - Send initialization commands to instrument. 

ADDR - GPIB address ( integer ) 

Does Selected Device Clear and Remote ENable. Should 
be called for every instrument intended for use. 


SUBROUTINES FOR TEKTRONIX 7912 DIGITIZER 


GETBESTMAI ( ADDR, MAI ) - Get the optimum intensity value 

MAI - Returns intensity value ( integer ) 

The number of interpolated points is found at current 
intensity; if 0, intensity is increased. 

Intensity is lowered step by step until the number 
of interpolated points reaches 0. This should 
be called before receiving data from 7912. 


GETMAI ( ADDR, MAI ) - Get current value of intensity from 7912 
MAI - Returns current intensity value ( integer ) 


SETVOLTDIV ( ADDR, SETTING ) - Set volts/div on 7912 

SETTING - Must be one of: .01, .02, .05, .1, .2, .5, 1., 2., 5. 


SETTIMEDI V( ADDR, SETTING ) - Set seconds/div on 7912 

SETTING - Must be one of: 1., 2., or 5. 

E -3, -6, -9, or -12 ( e.g. 2.E-6 ) 


AUTOS ETVOLTS ( ADDR, MAI, SETTING ) - Find best voltage/div setting 

MAI - Intensity to use ( integer ) 

SETTING - Returns volts/div setting ( real ) 

The lowest voltage setting ( .01 ) is tried initially 
and number of interpolated points on the 7912 is found. 

The setting is increased until interpolated points is 0. 

MAI should be value returned from GETBESTMAI. 


DIGDEF( ADDR ) - Do a "Digitize Defects" 

This operation should be done before taking data. It 
causes the 7912 to remember the defective points and 
does not send them when data is requested. 



GETGRID( ADDR, TIME, VOLTS ) - Get Time/division and Volts/Division 

TIME - Returns Seconds per division as set on the 7912 ( real ) 
VOLTS - Returns Volts per division as set on the 7912 ( real ) 


GETSA( ADDR, NAVE, MAI, A ) - Acquire Sample-Averaged waveform 

NAVE - Number of averages desired ( integer; 0 thru 64 ) 

MAI - Intensity to use ( integer ) 

A(512) - Integer array returned of added samples 

This subroutine sets intensity ( should be value 
returned from GETBESTMAI ), requests sample averaging, 
and requests read. If an error ocurrs, TEKRESET 
is called and read is requested again. The array A must 
be divided by NAVE to result in ( 1-255 ) - bottom to 
top of oscilloscope display. 


TEKGTL ( ADDR ) - Send GO TO LOCAL and TV MODE commands to 7912 

TEKRESET ( ADDR ) - Reset 7912 

This calls INITINSTR, GETBESTMAI, and DIGDEF 


SUBROUTINES FOR HEWLETT-PACKARD TIME SYNTHESIZER 


PUTTIME ( ADDR, TIME ) - Set Time Delay 

TIME - Desired Delay Time in seconds ( 0. < TIME < 166. E-3 ) 


GETTIME ( ADDR, TIME ) - Get Delay 

TIME - Returns time delay currently displayed ( real ) 


GETSTATE ( ADDR, STATE ) - Get State of Time Synthesizer 

STATE ( 66 ) - Byte array of codes describing the current 
complete state of the device. Useful for 
returning it back to a particular state with 
PUTS TATE. 


PUTSTATE ( ADDR, STATE ) - Set State of Time Synthesizer 

STATE ( 66 ) - A byte array of codes describing a state 
to set the device to. { See GETSTATE ) 


li 


SUBROUTINES FOR KLINGER XYZ CONTROLLER 


SETXYZ { ADDRX , IX, ADDRY, IY, ADDRZ , IZ ) - Initialize controller 

Before issuing any commands to the Klinger, this subroutine 
must be called to indicate the stepsize needed for each axis. 

IX - X axis stepsize, in microns ( integer; 0 < IX < 10000 ) 
IY - Y axis stepsize " IY 

IZ - Z axis stepsize " IZ 


MOVXYZ ( ADDR , IAXIS, IPM ) - Move one step 

I AX IS - Axis to move the distance according to most 
recent call of SETXYZ ( integer; [ 1,2,3 ]-[ X,Y,Z ] ). 

IPM - Flag for direction ( integer; <0 -backward, >0 -forward ) 


SUBROUTINES FOR FLUKE MULTIMETER 


GETFLUKE ( ADDR, VOLTAGE ) - Get voltage from Fluke Multimeter 
VOLTAGE - Value read by meter in volts. ( real ) 


SUBROUTINES FOR TEKTRONIX 7854 DIGITAL OSCILLOSCOPE 


GETWFM( ADDR, AVES, DATA, XINCR, YZERO, YSCALE ) - Acquire 

512-point averaged currently displayed waveform 

AVES - Number of averages ( integer ) 

DATA (512) - Real array representing averaged waveform 
XINCR - Time increment in seconds between points ( real ) 

YZERO - Value in volts represented by DATA value of 0 ( real ) 
YSCALE - Scale factor of DATA ( real; VOLTS/DATA value ) 


TEK78RESET( ADDR ) - Reset 7854 

This calls INITINSTR and does a serial poll. 


SUBROUTINES FOR TESTECH LS-86 SCANNER 


MOVTESTECH ( ADDR, IAXIS, IDEST ) - Move to destination IDEST 

IAXIS - Desired axis to move ( 1,2, 3, 4 )=( X,Y,Z, Swivel ) 
IDEST - Destination in thousandths of an inch if X,Y, or Z; 
hundredths of a degree if Swivel. ( integer ) 

X limits: 0 < IDEST < 36000 

Y limits: 0 < IDEST < 12000 

Z limits: 0 < IDEST < 15000 

Swivel limits: 0 < IDEST < 36000 


FINDTESTECH ( ADDR, 

IX, IY, 

IZ 

) - Get current location 

IX 

- Returns 

X 

location 

in 

thousandths ( 

integer ) 

IY 

- Returns 

Y 

location 

in 

thousandths 

ll 

IZ 

- Returns 

Z 

location 

in 

thousandths 

ft 


SUBROUTINES FOR HEWLETT-PACKARD A/D CONVERTER 


GETA2D ( ADDR, IDATA ) - Get value from channel 4 of H.P. A/D converter 

IDATA - Returns value from 0 to 1024, scaled such that 
1024 = approx 10 volts. ( integer ) 


13 


ononnn i-* h* noon 


Appendix C 

Low Level Source Code for Subroutines Listed in Appendix B 

c################################### 

c 

C GPIBASE . FOR 

C 

C This is a set of subroutines which communicate 

C with the IEEE-488 instruments found in the Structural 

C Integrity Branch laboratories. For further 

C documentation, see the file GPIBASE . TXT 

C 

C David Stang and Edward Generazio 

C 

C Latest update: 3-NOV-1988 

C 

C################################### 

SUBROUTINE STRTGPIB ( ICHAN ) 

Assign GPIB channel ( IEX unit ) 

Call IESTRT to initialize the unit 
Send IFC command to clear the interface 


byte 

integer*2 

integer*4 

character 

common 


WRKO ( 80 ) , BUFFER( 32 ) 
I STAT ( 4 ) 

ICHAN 

STRING*6 , I CHAR* 1 
/IBLK/ WRKO, BUFFER 


0020 


0100 


format( II ) 

ENCODE ( 1,10020, ICHAR ) ICHAN 
STRING = '_IXA'//ICHAR//' : ' 

ISTATUS - SYS$ASSIGN ( STRING , ICH ,,) 
WRITE ( 5,10100 ) ICH , ISTATUS 
CALL IESTRT( WRKO, ICH, 0,, 0, 1, ISTAT 
CALL IECMD( WRKO, ,, %REF('IFC'), ISTAT 


) 


FORMAT ( 

RETURN 

END 


CHANNEL ASSIGNED 


,15, 


STATUS 


=',I2 ) 


SUBROUTINE INITINSTR( LISTEN ) 

Call IEPOLL, send SDC, and any needed 
initialization commands to device 
at address LISTEN 


byte WRKO (80), BUFFER( 32 ) 
byte ADDR ( 2 ) , CHR8 

common /IBLK/ WRKO, BUFFER 
integer*2 ISTAT(4) 

NADDR =* 1 

ADDR ( 1 ) - LISTEN + 32 
ADDR ( 2 ) - 96 

C I F ( LISTEN. EQ.O )NADDR=2 

CALL IECMD( WRKO, ADDR, NADDR, 
CALL IECMD( WRKO, ADDR, NADDR, 
RETURN 
end 


! IF THERE IS A SECONDARY ADDRESS 
! TWO BYTES MUST BE SENT 
%REF ( ' SDC ' ) , ISTAT ) 

%REF ( ' REN ' ) , ISTAT ) 


C 
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SUBROUTINE GETSA( ADDR , NAVE, MAI, A ) 

Acquire 512-point waveform from Tek 7912 

NAVE - Number of averages 

MAI = Main Intensity 

A( 512 ) = Data 


byte 
byte 
character 
integer *2 
integer*4 
equivalence 
equivalence 
common 


BUFFER( 32 ) , WRKO(80) 

ABYTE ( 1029 ) , LADDR( 2 ) , TADDR( 2 ) 
AINT* 3 , ANAVE * 2 , ABUFFER* 32 
A( 512 ) , AWORD( 512 ) , ISTAT(4) 
ADDR 

( ABUFFER, BUFFER ) 

( AWORD(l), ABYTE ( 4 ) ) 

/IBLK/ WRKO, BUFFER 


B 


100 LADDR(l) = ADDR+32 
LADDR( 2 ) = 96 
TADDR(l) = ADDR+64 
TADDR ( 2 ) =96 


Set intensity 

ENCODE ( 3,10023, AINT )MAI 
3 format( 13 ) 

ABUFFER (1:8) = 'MAI '//AINT//';' 

CALL IESEND( WRKO, BUFFER, 8, LADDR , 2 , ISTAT ) 


Request Sample Averaging 

ENCODE ( 2,10022, ANAVE ) NAVE 
0022 format( 12 ) 

ABUFFER( 1:10) = 'DIG SA, '//ANAVE//' ; ' 

CALL IESEND( WRKO, BUFFER, 10, LADDR, 2, ISTAT ) 

Request Read, Get data, check for error ( e.g. timeout error ) 

ABUFFER (1:6) = 'REA SA' 

BUFFER( 7 ) = 13 

CALL IESEND( WRKO, BUFFER, 7, LADDR, 2, ISTAT ) 

CALL WAIT2{ 400) 

CALL IERECV { WRKO, ABYTE, 1029, ,, TADDR, ISTAT ) 

I F ( I STAT ( 1 ) . NE . 1 ) THEN 
WRITE ( 8,10800 ) 

10800 format( ' ERROR AT IERECV IN GETSA, TRYING AGAIN' ) 

CALL TEKRESET ( AD, MAI ) 

GOTO 100 
END IF 


DO 1000 1=1,512 ! 

II = 2*1 + 2 ! Byte swap ! 

B= ABYTE ( II ) 

ABYTE ( II ) = ABYTE { II+l ) 

ABYTE ( II+l ) = B 
1000 A( I ) = AWORD ( I ) 

RETURN 

end 

C 
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SUBROUTINE AUTOSETVOLTS ( AD, MAI, SET ) 

C 

C Sets appropriate voltage setting by 

C checking on number of interpolated points 

C 

C MAI - Intensity 

C SET - Returns setting 

C 

byte BUFFER( 32 ) , WRK0(80) 

character ABUFFER* 3 2 , AMAI*3 
byte ADDR ( 2 ) 

real SETTINGS(6) 

integer *2 ISTAT<4) 
integerM AD 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 

data SETTINGS/ 0.01,0.02,0.05,0.10,0.20,0.50 / 

ADDR ( 1 ) - AD+32 

ADDR( 2 ) = 96 

ENCODE ( 3,10023, AMAI )MAI 

10023 format( 13 ) 

1 = 1 
IB=0 

C TYPE *, 'CALLING SETVOLTSDIV' 

100 CALL SETVOLTDI V ( AD , SETTINGS ( I ) ) 

WRITE ( 5, 1133 )SETTINGS( I ) 

1133 FORMAT( '+', 'CURRENT GAIN',F10.5) 

CALL WAIT2( 100) 

C TYPE * , ' SETVOLTS ANSWERED ', SETTINGS ( I ) 

ABUFFER (1:26) = 'MAI '//AMAI//' ; DIG DAT ; DEF ON;ATC' 
CALL I ESEND ( WRKO, BUFFER, 26, ADDR, 2, ISTAT ) 

CALL WAIT2( 100 ) 

CALL GETNUMINT ( AD , NUMINT , MAI ) 

C TYPE *, NUMINT 

I F ( NUMINT. EQ. 999 ) THEN 
SET = 999. 

GOTO 900 
ENDIF 

IF ( NUMINT. EQ.O ) THEN 

I F ( IB.EQ.l ) GOTO 200 
IB=1 

GOTO 100 
ENDIF 

200 I F ( NUMINT. GT . 0 ) THEN 

1 = 1 + 1 
GOTO 100 
ENDIF 

SET = SETTINGS (I) 

900 RETURN 

end 
C 

SUBROUTINE GETBESTMAI ( AD, MAI ) 

Figure out the best intensity by finding number 
of interpolated points over a number of trials 

byte BUFFER( 32 ) , WRK0(80) 

character ABUFFER* 32 , AMAI*3 
byte ADDR( 2 ) 

intege r *2 ISTAT(4) 
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10023 

C 

C 

C 

100 

C 


900 

C 
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C 

10023 


integer*4 AD, MAI 
equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 
format( 13 ) 

ADDR(l) = AD+32 
ADDR ( 2 ) = 96 
IB - 0 

TYPE * , ' GONE TO LOCAL IN GBM ' 

CALL GETMAI( AD , I TRY ) 

TYPE * , ' GOT MAI AGAIN IN GET BESTMAI ' 

CALL DIGDEF ( AD ) 

CALL WAIT2 ( 2000 ) 

ENCODE ( 3,10023, AMAI ) I TRY 

ABUFFER( 1:26) = 'MAI '//AMAI//' ;DIG DAT ; DEF ON;ATC' 
CALL IESEND( WRKO, BUFFER, 26, ADDR, 2, ISTAT ) 

CALL WAIT2 ( 100 ) 

CALL GETNUMINT ( AD, NUMINT, MAI ) 

TYPE *, NUMINT 
I F ( NUMINT. GT.O ) THEN 
ITRY-ITRY + 1 
IB - 1 
GOTO 100 
ENDIF 

I F ( IB.EQ.l ) GOTO 900 
ITRY=ITRY - 1 
IB - 0 
GOTO 100 
MAI = ITRY + 20 
RETURN 
end 


SUBROUTINE GETMAI ( AD, MAI ) 

Get intensity from TEK 7912 

byte BUFFER( 32 ) , WRKO (80) 

character ABUFFER*32 
byte ADDR ( 2 ) 

integer*2 ISTAT(4) 

integer*4 AD, MAI 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 
ADDR ( 1 ) - AD+32 
ADDR ( 2 ) - 96 
ABUFFER( 1:4) - 'MAI?' 

CALL IESEND ( WRKO, BUFFER, 4, ADDR, 2, ISTAT ) 
ADDR ( 1 ) - 64 

CALL IERECV( WRKO, BUFFER, 20, ,, ADDR, ISTAT ) 
OUT-O 

DO 5 K-5,20 

I F ( OUT . EQ . 1 ) GOTO 5 

I F ( BUFFER(K) . EQ . ' ; ' ) THEN 

KK-K 

OUT-1 

ENDIF 

CONTINUE 

DECODE ( KK-5, 10023, ABUFFER ( 5 : KK-1 )) MAI 
DECODE ( 20,10023, ABUFFER( 5:20) )MAI 
format( 13 ) 

WRITE ( 5,1133 )MAI 


non 


1133 


10020 

40 

100 

10100 

200 

10200 

300 

1133 

C 


FORMAT ( ' + ' , 'CURRENT INTENSITY' , 17 ) 

RETURN 

end 

SUBROUTINE GETNUMINT ( AD, NUM, MAI ) 

Get number of interpolated points 

byte BUFFER( 32 ) , WRK0(80) 

character ABUFFER*32 
byte ADDR( 2 ) 

integer*2 ISTAT(4) 

integer*4 AD, NUM, MAI 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRK0, BUFFER 
format( I ) 

ADDR(l) - AD+32 
ADDR( 2 ) - 96 
IB-0 

ABUFFER( 1:4) = 'INT?' 

CALL IESEND( WRK0, BUFFER, 4, ADDR, 2 , ISTAT ) 

ADDR(l) - 64 
CALL WAIT2 ( 100 ) 

CALL IERECV( WRK0, BUFFER, 20, ,, ADDR, ISTAT ) 

I F ( BUFFER ( 5 ) . EQ . ' N ' ) THEN 

I F ( IB.GT.O ) THEN 

WRITE ( 8,10100 ) 

format( ' "N" RECIEVED IN GETNUMINT, RESETTING') 

CALL TEKRESET( AD , MAI ) 

GOTO 40 
ENDIF 

IB-1 

GOTO 100 
ENDIF 

K-l 

I F ( BUFFER(K) .EQ. '?' ) THEN 

WRITE ( 8,10200 ) 

format ( ' "?" RECIEVED IN GETNUMINT, RESETTING') 

CALL TEKRESET ( AD, MAI ) 

GOTO 40 

ENDIF 

I F ( BUFFER( K) .EQ. ' ; ' ) THEN 

DECODE ( K-4, 10020, ABUFFER( 4 : K-l ) )NUM 
GOTO 300 
ENDIF 
K-K+l 
GOTO 200 

WRITE (5,1133 )NUM 

FORMAT ('+',' INTERPOLATED POINTS ',17) 

RETURN 

end 


SUBROUTINE SETVOLTDIV( 


AD, SETTING ) 


byte 

character 

character 

byte 

integer*2 


BUFFER ( 32 ) , 
ABUFFER*32 
AVOLTS*10 
ADDR ( 2 ) 

ISTAT( 4 ) 


WRK0 (80 ) 
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integer*4 AD 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 
ADDR(l) - AD+32 
ADDR ( 2 ) = 97 
10040 format( E10.5 ) 

ENCODE ( 10, 10040, AVOLTS ) SETTING 

ABUFFER (1:14) = 'V/D '//AVOLTS 

CALL I ESEND ( WRKO, BUFFER, 14, ADDR, 2, ISTAT ) 

RETURN 

end 

SUBROUTINE SETTIMEDIV ( AD, SETTING ) 

byte BUFFER (32), WRK0(80) 

character ABUFFER*32 
character ATIME*10 
byte ADDR ( 2 ) 

integer *2 ISTAT(4) 
integer*4 AD 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 

ADDR(l) - AD+32 
ADDR ( 2 ) = 98 
10040 format( E10.5 ) 

ENCODE ( 10, 10040, ATIME ) SETTING 
ABUFFER( 1:14) - 'T/D '//ATIME 

CALL I ESEND ( WRKO, BUFFER, 14, ADDR, 2, ISTAT ) 

RETURN 

end 

SUBROUTINE DIGDEF ( AD ) 

byte BUFFER( 32 ) , WRK0(80) 

character ABUFFER* 32 
byte ADDR ( 2 ) 

integer*4 ad 
integer*2 ISTAT(4) 
equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 

ADDR ( 1 ) - AD+32 
ADDR ( 2 ) - 96 

ABUFFER( 1:11) *= 'DIG DEF,100' 

CALL I ESEND ( WRKO, BUFFER, 11, ADDR, 2, ISTAT ) 
CALL WAIT2 ( 2000 ) 

RETURN 

end 


SUBROUTINE GETGRID( AD , TIME , VOLTS ) 

Get Time/Div and Volts/Div from TEK 7912 

byte BUFFER (32), WRK0(80) 

character ABUFFER* 32 
byte ADDR( 2 ) 

integer*2 ISTAT(4) 

integer*4 AD 

equivalence ( BUFFER, ABUFFER ) 
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10040 


200 


300 

400 

900 

C 


common /IBLK/ WRKO, BUFFER 

ADDR(l) = AD+32 
ADDR ( 2 ) = 96 
format( E ) 

ABUFFER( 1:4) = 'HSl?' 

CALL IESEND( WRKO, BUFFER, 4, ADDR, 2, ISTAT ) 

ADDR(l) =64 

K=1 

CALL IERECV ( WRKO, BUFFER, 15, ,, ADDR, ISTAT ) 

I F ( BUFFER(K) . EQ . ' ; ' ) THEN 

DECODE ( K-4, 10040, ABUFFER( 4 : K— 1 ) ) TIME 
GOTO 300 
END IF 
K=K+1 
GOTO 200 

ADDR(l) = 32 
ABUFFER (1:4) = 'VSl?' 

CALL IESEND( WRKO, BUFFER, 4, ADDR, 2, ISTAT ) 

ADDR ( 1 ) = 64 

K=1 

CALL IERECV( WRKO, BUFFER, 15, ,, ADDR, ISTAT ) 

I F ( BUFFER(K) .EQ. ' ; ' ) THEN 
DECODE ( K-4, 10040, ABUFFER( 4 : K— 1 ) )VOLTS 
GOTO 900 
ENDIF 
K=K+1 
GOTO 400 
RETURN 
end 


SUBROUTINE TEKGTL ( AD ) 

Tell 7912 to Go To Local & put into TV MODE 

byte BUFFER( 32 ) , WRKO (80) 

character ABUFFER* 32 
byte ADDR ( 2 ) 

integer*2 ISTAT(4) 

integer*4 AD 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 
ADDR(l) = AD+32 
ADDR( 2 ) = 96 

ABUFFER( 1:8) = 'MODE TV;' 

CALL IESEND( WRKO, BUFFER, 8, ADDR, 2, ISTAT ) 
CALL IECMD( WRKO, ADDR, 2, IREF('GTL'), ISTAT ) 
RETURN 
end 

SUBROUTINE TEKRESET ( AD, MAI ) 

byte BUFFER( 32 ) , WRK0(80) 
integer*4 AD, MAI 
common /IBLK/ WRKO, BUFFER 
CALL INITINSTR ( 0 ) 

CALL SETVOLTDIV ( AD,. 5 ) 

CALL GETMAI( AD, MAI ) 

CALL GETBESTMAI ( AD, MAI ) 
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c 


CALL DIGDEF ( AD ) 

RETURN 

end 


C 


SUBROUTINE PUTTIME( AD, TIME ) 

C 

C Set time synthesizer delay to TIME 

C 

byte BUFFER(32), WRK0{80) 
character CTIME*10 
integer*2 ISTAT{4) 
integer*4 AD 
byte ADDR 

equivalence ( BUFFER ( 2 ) , CTIME ) 
common /IBLK/ WRKO, BUFFER 

ADDR - AD+32 

ENCODE ( 10,10100, CTIME ) TIME 
10100 format( E10.3 ) 

BUFFER ( 1 ) - 'D' 

BUFFER( 12 ) = 

CALL IESEND( WRKO, BUFFER, 12, ADDR , 1 , ISTAT ) 

RETURN 

end 


SUBROUTINE GETTIME ( AD , D ) 

C 

C Get Time delay from synth 

C 

byte BUFFER( 32 ) , WRKO (80), STATE(66) 
byte ADDR 

integer*2 DAVE(66), ISTAT(4) 
integer*4 AD 

real U(12), P ( 12 ) , G(7) 

common /IBLK/ WRKO, BUFFER 
CALL GETDAVE ( AD, DAVE ) 

U( 12 ) - REAL(DAVE(12) ) 

B = 0. 

DO 200 1= 6,11 
U(I) - REAL ( DAVE ( I ) ) 

P(I) - U(I) - 6 . *REAL { ( I I NT ( U( I )/16 . ) 
G( 1-5 ) =* P(I)*( 10 . ** ( -2* ( 1-5 ) )) 

200 B=B + G( 1-5 ) 

A - U( 12 )-244 . 

D - B * ( 10.**( A-12 . ) ) 

RETURN 
end 



The algorithm for converting the 
array DAVE(66) to the time delay D 
was developed by Alex Vary, 
NASA Lewis Research Center, 
Cleveland, Ohio. 


c 

SUBROUTINE PUTSTATE ( AD, STATE ) 

C 

C Put State into Time Synth 

C 

byte BUFFER( 32 ) , WRKO (80), STATE ( 66 ), ADDR 
integer*2 ISTAT(4) 
integer*4 AD 

common /IBLK/ WRKO, BUFFER 
ADDR - AD+32 
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BUFFER ( 1 ) - 'D' 

BUFFER(2) - 13 

CALL IESEND( WRKO , BUFFER, 2, ADDR , 1 , ISTAT ) 
BUFFER(l) - 'L' 

BUFFER( 2 ) - 'N' 

CALL IESEND( WRKO, BUFFER, 2, AODR , 1 , ISTAT ) 
CALL IESEND( WRKO, STATE, 66, ADDR, 1 , ISTAT ) 
RETURN 
end 

SUBROUTINE GETSTATE ( AD, STATE ) 

Get State from Time Synth 

byte BUFFER( 32 ) , WRKO (80), STATE(66) 
byte LADDR , TAD DR 

integer*2 ISTAT(4) 
integer*4 AD 

common /IBLK/ WRKO, BUFFER 

TADDR - AD+32 
LADDR - AD+64 
BUFFER(l) - 'D' 

BUFFER( 2 ) - 13 

CALL IESEND( WRKO, BUFFER, 2, TADDR, 1, ISTAT ) 
BUFFER ( 1 ) - 'T' 

BUFFER( 2 ) - 'E' 

BUFFER ( 3 ) - 13 

CALL IESEND( WRKO, BUFFER, 3, TADDR, 1, ISTAT ) 
CALL I ERECV ( WRKO, STATE, 66, ,, LADDR, ISTAT ) 

RETURN 

end 


SUBROUTINE GETDAVE ( AD, DAVE ) 

Get State from Time Synth as integers with stuff 
in the MSB 

byte BUFFER( 32 ) , WRKO (80), STATE (66), B 
byte TADDR, LADDR 
integer*2 DAVE(66), ISTAT(4), I 
integer*4 ad 
character MSB*2, LSB*1 
equivalence ( MSB, I ) 
equivalence ( LSB , B ) 
common /IBLK/ WRKO, BUFFER 

TADDR - AD+32 
LADDR - AD+64 
BUFFER(l) - 'D' 

BUFFER( 2 ) = 13 

CALL IESEND ( WRKO, BUFFER, 2, TADDR, 1, ISTAT ) 
BUFFER(l) = 'T' 

BUFFER( 2 ) - 'E' 

BUFFER ( 3 ) = 13 

CALL IESEND ( WRKO, BUFFER, 3, TADDR, 1, ISTAT ) 

CALL I ERECV ( WRKO, STATE, 66, ,, LADDR, ISTAT ) 

1-0 

K=0 
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DO 700 J-1,20 

B = STATE (J+K) 

MSB (1:1)= LSB 
DAVE (J+K) = I 
B = STATE (J+K+l) 

MSB (1:1)= LSB 
DAVE (J+K+l) = I 
700 K=K + 1 
RETURN 
end 


SUBROUTINE SETXYZ( AD,Il, A2 , 12, A3, 1 3 ) 

Initialize Klinger with step values 

byte BUFFER( 32 ) , WRK0(80) 
byte ADDR , BSTEP ( 6 ) 
integer ISTEP( 3 ) , AD , Al , A2 , A3 

integer*2 ISTAT(4) 
character ABUFFER*25, ASTEP*6 
equivalence ( ABUFFER, BUFFER ) 
equivalence ( ASTEP, BSTEP ) 
common /IBLK/ WRKO, BUFFER 

format( I ) 
format( 16 ) 
format( ' ' ,A32 ) 

ISTEP(l) - II 
I STEP ( 2 ) - 12 
ISTEP ( 3 ) - 13 
ABUFFER( 1:5) - 'R 253' 

BUFFER ( 6 ) =13 

ABUFFER (7:9) = 'S 2' 

BUFFER(IO) =13 
ABUFFER ( 11:14) ='F 20' 

BUFFER (15) =13 
BUFFER( 16 ) = 'N' 

BUFFER( 17 ) = ' ' 

DO 200 1=1,3 

NDIG = I FIX ( ALOG1 0 ( REAL ( ISTEP ( I ) ) ) )+l 

N6M » 6-NDIG 

ENCODE ( 6,10026, ASTEP )ISTEP(I) 

DO 90 IDIG=1 , NDIG 
90 BUFFER ( 17+IDIG) = BSTEP ( N6M+IDIG ) 

BUFFER ( NDIG+18 ) = 13 
BUFFER ( NDIG+19 ) = 'A' 

BUFFER(NDIG+20 ) = 13 
NSEND = 20+NDIG 
ADDR = AD+32+( 1-1 ) 

200 CALL IESEND( WRKO, BUFFER, NSEND , ADDR, 1 , ISTAT ) 
RETURN 
end 

SUBROUTINE MOVXYZ ( AD,IAXIS, I PORN ) 

C 

C Move Klinger one step 

C 

byte BUFFER( 32 ) , WRKO (80), ADDR 
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C 

C 

c 


10020 

10026 

10124 


Here is where we stick STATE 
into the MOST SIGNIFICANT BYTE 
of DAVE 
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integer*2 ISTAT(4), LISTEN 
integer*4 AD 

common /IBLK/ WRKO, BUFFER 

ADDR - AD + 32 

I F ( IPORN.GT.O ) BUFFER( 1 ) = '+' 

I F ( IPORN.LT. 0 ) BUFFER ( 1 ) = 

BUFFER( 2 ) = 13 
BUFFER ( 3 ) = 'G' 

BUFFER( 4 ) = 13 

CALL I ESEND ( WRKO, BUFFER, 4, ADDR , 1 , ISTAT ) 

RETURN 

end 

C 

SUBROUTINE GETFLUKE ( AD, VOLTS ) 

Get value from Fluke Multimeter 

byte BUFFER( 32 ) , WRKO (80) 

character ABUFFER*32 
byte ADDR 

integer*2 ISTAT(4) 

integer*4 AD 

equivalence ( BUFFER, ABUFFER ) 
common /IBLK/ WRKO, BUFFER 
10 ADDR = AD+32 

ABUFFER (1:3) = 'VC?' 

CALL IESEND( WRKO, BUFFER, 3, ADDR, 1 , ISTAT ) 

ADDR - 69 

CALL IERECV( WRKO, BUFFER, 14, ,, ADDR, ISTAT ) 

I F ( BUFFER(l) .NE. '+' .AND. BUFFER ( 1 ). NE ' )GOTO 10 
DECODE ( 12,10040, ABUFFER( 1:12) )VOLTS 
0040 format( E ) 

RETURN 

end 

SUBROUTINE WAIT2(NWAIT) 

DO 1000 IWAIT=1 , NWAIT 
DO 500 1=1, NWAIT 
do 501 j=l,nwait 
501 continue 
500 CONTINUE 
1000 CONTINUE 
RETURN 
END 


SUBROUTINE MOVTESTECH ( AD, IXYORZ, IDEST ) 

Move TESTECH scanner IXYORZ axis ( 1,2, 3, 4 )=( X,Y,Z,S ) 
to destination IDEST 

BYTE BUFFER( 48 ) , WRKO (80) 

BYTE IXYORZ , B , B6 ( 6 ) ,SCODE(4) , CR , TADDR , LADDR 

INTEGER* 4 AD 
CHARACTER* 6 C6 
EQUIVALENCE ( B6,C6 ) 

COMMON /IBLK/ WRKO, BUFFER 

DATA SCODE/ 68,70,87,73 /, CR/13/ 
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TAD DR - AD+32 
LADDR - AD+32 
I F ( IDEST.LT. 0 ) THEN 

TYPE *,' ILLEGAL DESTINATION ' 

GOTO 900 
ENDIF 

BUFFER(l) = SCODE ( IXYORZ ) 

Convert IDEST to a 6-byte character string 
with leading zeros and place in BUFFER(2:7) 

ENCODE ( 6,10070, C6 ) IDEST 
DO 100 1=1,6 
B=B6 ( I ) 

I F ( B.EQ.32 ) B=48 
100 BUFFER( 1+1 ) = B 

1070 FORMAT ( 16 ) 

BUFFER( 8 ) = CR 

Call IESEND to move 

CALL IESEND( WRK0, BUFFER, 8, TADDR , 1 , ISTAT ) 

CALL IERECV ( WRK0, BUFFER, 2, ,, LADDR, ISTAT ) 

I F ( IXYORZ. EQ. 4 ) THEN 
BUFFER ( 1 ) = 'Z' 

BUFFER( 2 ) = CR 

CALL IESEND( WRK0, BUFFER, 2, TADDR, 1, ISTAT ) 

CALL IERECV ( WRK0, BUFFER, 2, ,, LADDR, 1, ISTAT ) 

ENDIF 

900 RETURN 
END 

SUBROUTINE FINDTESTECH ( AD , IXPOS , IYPOS , IZPOS ) 

Ask TESTECH scanner for current position and 
return X, Y, and z axis positions as integers 
in units of thousandhs of an inch 

BYTE WRKl ( 80 ) , BUFFER<48) 

BYTE GETPOS ( 62 ) , CR 
INTEGER*2 ISTAT(4) 

INTEGER* 4 AD 

CHARACTER*6 XPOS , YPOS , ZPOS 

EQUIVALENCE ( XPOS, GETPOS ( 1 ) ), ( YPOS, GETPOS ( 7 ) ) 

EQUIVALENCE ( ZPOS, GETPOS (13) ) 

COMMON /IBLK/ WRKl, BUFFER 
DATA CR/13/ 

TADDR = AD+32 
LADDR = AD+64 
BUFFER(l) = ' K ' 

BUFFER( 2 ) = CR 

CALL I ESEND ( WRKl, BUFFER, 2, TADDR, 1, ISTAT ) 

CALL I ERECV ( WRKl, GETPOS, 62, ,, LADDR, ISTAT ) 

I F ( ISTAT(l) .NE.l ) CALL LIB$SIGNAL( %VAL( ISTAT( 1 ) ) ) 

DECODE ( 6,10060, XPOS ) IXPOS 
DECODE ( 6,10060, YPOS ) IYPOS 
DECODE ( 6,10060, ZPOS )IZPOS 
10060 FORMAT ( 16 ) 
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RETURN 

END 


SUBROUTINE GETA2D( AD, DATA ) 

Get one word of data from channel 4 of the H.P. A/D converter 

BYTE WRKO ( 80 ) , BUFFER( 48 ) 

BYTE A( 2 ) , B ( 2 ) , TADDR , LADDR 

INTEGER*2 D , DATA , ISTAT(4) 

INTEGER* 4 AD 
EQUIVALENCE ( B , D ) 

COMMON /IBLK/ WRKO, BUFFER 
0010 FORMAT ( A3 ) 

0020 FORMAT ( A2 ) 

TADDR = AD+32 

LADDR = AD+64 

Send "H8A", wait, then send "JF" 

ENCODE ( 6,10010, BUFFER )'H8A' 

CALL IESEND( WRKO, BUFFER, 3, TADDR, 1, ISTAT ) 

CALL WAIT 3 ( 1 ) 

ENCODE ( 6,10020, BUFFER )'JF' 

CALL IESEND ( WRKO, BUFFER, 2, TADDR, 1, ISTAT ) 

Get data, put 1st byte in DATA'S 2nd byte, etc 

CALL I ERECV ( WRKO, A, 2, ,, LADDR, ISTAT ) 

B ( 1 ) - A(2) 

B(2) = A( 1 ) 

DATA - D 

RETURN 

END 

SUBROUTINE WAIT3(NWAIT) 

DO 1000 IWAIT=1 , NWAIT 
DO 500 1=1,1000 
500 CONTINUE 
1000 CONTINUE 
RETURN 
END 
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